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Abstract 

The values at pseudothreshold of two loop sunrise master amplitudes with arbi- 
trary masses are obtained by solving a system of differential equations. The expan- 
sion at pseudothreshold of the amplitudes is constructed and some lowest terms are 
explicitly presented. 



PACS ll.lO.-z Field theory 

PACS ll.lO.Kk Field theories in dimensions other than four 
PACS ll.15.Bt General properties of perturbation theory 



1 Introduction. 



The sunrise graph (also known as sunset or London transport diagram) appears naturally, 
as a consequence of tensorial reduction, in several higher order calculations in gauge 
theories. Due to the presence of heavy quarks, vector bosons and Higgs particles all the 
internal lines may carry a different mass, so that sunrise amplitudes depend in the general 
case on three internal masses rrii, i = 1,2, 3, and the external scalar variable p^, if is 
the external momentum (in n-dimensional Euclidean space). 

For a proper understanding of their behaviour as well as for a check of the numerical 
calculations, it is convenient to know the amplitudes off-shell, and also around some par- 
ticular values of p^, such as p"^ = 0, p"^ = oo, at the threshold and at the pseudothresholds. 

This paper is devoted to the analytic evaluation of the coefficients of the expansion of 
the sunrise amplitudes in p"^, at the pseudothreshold value p"^ = —{mi + m2 — 
( the other pseudothreshold values —{mi — m2 + m^Y, —{mi — m2 — m^Y, can be easily 
obtained by permutation of the masses ); the approach rehes on the exploitation of the 
information contained in the linear system of first order differential equations in p^, which 
is known to be satisfied by the sunrise amplitudes themselves It is to be noted that 
all the above points ( = 0, oo , threshold and pseudothresholds) correspond to the 
Fuchsian points of the differential equations, which therefore emerge as a natural tool for 
their discussion. 

The analytic properties of Feynman diagrams at threshold and pseudothresholds are 
well known, see for example . The sunrise diagram, with different masses, has been in- 
vestigated in while in [Q] the values of the amplitudes at threshold and pseudothreshold 
were obtained. With the method established in Q and , further, the expansion around 
threshold was recently obtained in 0. 

The sunrise amplitudes are regular at any of the pseudothresholds, say Sq, so that 
they can be expanded as a single power series in {p^ — sq) around that point [§]. When 
the expansions are inserted in the differential equations, the equations become a set of 
algebraic equations in the coefficients of the expansions; the obtained algebraic equations 
can then be solved recursively, for arbitrary value of the dimension n, once the initial 
conditions, i.e. the values of the sunrise amplitudes at the considered pseudothreshold, 
are given. As those initial values are in turn functions of the masses, we find in our 
approach a system of linear differential equations in the masses satisfied by them. We 
expand the equations in the masses in the dimension n around n = 4 and solve them 
exphcitly up to the finite part in (n — 4). This result for n = 4 and p^ = Sq is in 
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agreement with the hterature [^. 

We then look at the recursive solution of the algebraic equations for the coefficients 
of the expansion in (p^ — Sq), not yet discussed in the literature for the pseudothreshold 
(our results are given in the Sections 4 and 5). It turns out that the formula, expressing 
the first order coefficients of the (p^ — Sq) expansion in terms of the zeroth order values, 
involves the coefficient 1/ (n — 4) . Therefore the finite part at {n — 4) of the first order 
terms in the {p^ — Sq) expansion involve the first order terms in (n — 4) of the zeroth 
order terms in (p^ — sq). Rather than evaluating those first order terms in (ra — 4), we 
prefer to evaluate directly the required finite part in n of the first order terms of the 
(p^ — So) expansion by means of a suitably subtracted dispersion relation in (p^ — sq) at 
n = 4. Fortunately, the formulae expressing the higher order coefficients of the {p^ — sq) 
expansion involve coefficients like l/(n — 5), 1/(?t, — 6) etc., which are finite at n = 4, and 
can be used without further problems for evaluating those higher order terms. 

The plan of the paper is as follows. In the second section the system of differential 
equations for the sunrise master amplitudes is recalled; in the third section the resulting 
differential equations in the masses for the values at the pseudothreshold are derived, 
expanded in [n — 4) and solved explicitly up to the finite part in (n — 4); in the fourth 
section the expansion in at the pseudothreshold is discussed; in the fifth section the 
finite part in {n — 4) of the first term in the expansion is given. The sixth section 
is a short summary, while the Appendix contains a list of various functions (mainly 
polynomials) used in the text. 



2 The equations. 

It is known that the two-loop sunrise self- mass graph with arbitrary masses mi, 7712,7713 
has four independent master amplitudes 0, which are referred to, as in |T|, by 

Fa{n,rnl,7nl,7nl,p^), a = 0,1,2,3, (1) 

where n is the continuous number of dimensions, mj,i = 1,2,3 the three masses and 
the external n-momentum. FQ{7i,7n\,7n\,7n\,p'^) is the scalar amplitude 

Fo{7i,7n^,7n^,7n^,p ) = 
J {2n)^-^J (27r)"-2 (^kf + 77il)m + 7nl){{p ~ h - + 771I) ' ^ ^ 



2 



while for z = 1, 2, 3 



771 / 2 2 2 2\ 



d 



771/ 2 2 2 2\ 



(3) 



but the 4 amphtudes are otherwise independent {i.e. none of them can be expressed as a 
hnear combinations of the others times rational functions of the arguments). 

In it was shown how to write for them a system of first order linear differential 
equations of the form 



d 



P^^^ol"-, ml, ml, p^)=(n - 3)Fo(n, ml, ml, ml,p'^) + ^ mlFi{n, ml, ml, mj, p 



i=l 



d 



2 f 2 2 2 2\ 

P -Q^Fi{n,m^,m2,m^,p 



EA / 2 2 2 2\77/ 2 2 2 2\iD/ 2 2 2 2\ 

Ai^a{n,m^,m2,m^,p )Fa{n,m^,m2,m^,p ) + Bi{n,m^,m2,m^,p ) 



a=0 



(4) 



where the Ai^^in, ml, ml, m\,p^) are known polynomials of the arguments, and the func- 
tions Bi{n,ml,ml,m1,p'^), constituting the non homogeneous part of the equations, are 
also known, being the combinations of other known polynomials times the products of two 
one-denominator, one-loop vacuum amplitudes of two different masses, T{n, m1)T{n, m^). 



where, as in 



T{n, m^) = C{n)- 



m 



n-2 



(n-2)(n-4) ' 

and C{n) is a coefficient depending on n which at n = 4 takes the value C(4) = 1 (its 
explicit form, which is essentially irrelevant in practice, can be found in [Q). 

The explicit form of the Aia{n,ml,ml,ml,p'^),Bi{n,ml,ml,ml,p'^), appearing in 
Eq.(^ is given in and is not repeated here for the sake of brevity; for the follow- 
ing, it is sufficient to give the explicit expression of the polynomial D{ml,ml,m\,p'^), 
which reads 



7^/ 2 2 2 2\ 

D{m^,m2,m^,p 



p^ + (mi + m2 + m-iY p^ + (mi + m2 - ms)^ 
p^ + (mi — m2 + ma)^ p^ + (mi — m2 — ms)^ 



(6) 



The above polynomial vanish when p^ takes one of the three pseudothreshold values 
— (mi -|- m2 — Tn^Y, —(mi — m2 + m^Y, —(mi — m2 — fn^Y, and at the physical threshold 
p^ = —{mi + m2 + msY- It is known that the Fj(n, mf, m^, m|, p^) and their p^-derivatives 



3 



are regular at the pseudothresholds; therefore, the apparent pseudothreshold pole in the 
r.h.s. of Eq.(^) must be canceled by a corresponding zero in the numerator. 



The recurrence relations solving the integration by part identities of |]T0[ for the sunrise 
amplitudes can be also written (see or the Appendix of [Q) in a form similar to Eq.(|), 
such as 



d 



^-^ I 2 2 2 2 



1 



ED / 2 2 2 2\77/ 2 2 2 2\ic/ 2 2 2 2\ 



a=0 



(7) 



where i = 1,2,3, the Ri^a{n,rn\^m\,m\,p'^) and the Si{n,m\,m\,m\^p'^) have the same 
structure as the y4j^Q,(n, , m|, m|,p^) and the Bi{n,m\,m\^m\,p'^) of Eq.(^; as in 
Eq.(H), the apparent pseudothreshold singularity of Eq.(|^) must be canceled by a cor- 
responding zero of the numerator. 



3 The pseudothreshold values of the master ampli- 
tudes. 

As the master amplitudes are regular at the pseudothresholds, considering for definiteness 
the pseudothreshold at p^ = —{mi + m2 — m^)'^, we can expand Fo(n, m^, m|, mg,^^) in 
p"^ around that point as 

Fa{n,ml,ml,ml,p'^) = if(°'°)(n, mi, mg, mg) 

+ H^"'^\n, mi, 7712, {j? + (777i + ?772 - "^3)^^ 

+ (77, 7771,7772,^3) (p^ + (mi + 7712 - "^3)^) + 

a = 0,1,2,3, (8) 

from which we can write at once, by using the definitions Eq.(|^), 

Go{n,mi,m2,m3) = Fo{n,m\,ml,ml,p^ = -{mi + m2 - m^Y) 

= if^°'°^(77, mi, m2, m3) , 

Gi(77, mi,m2,m3) = Fi{n,m\,m\,ml,p^ = -{mi + m2 - m^Y^ 

= H^^''^\n,mi,m2,m3) 

9 r7(0,0)/ N "^1 +"^2 -"^3 rr(o,l)/ N 

-U^ ' '[n,mi,m2,m^) ' '[n,mi,m2,m^) , 



dml mi 
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^(2,0), 



71,7711,1712,1713) 



9 rj(0,0)/ N "^1 +"^2 -"^3 r7(0,l)/ N 
— -U^ ' '[71,7171, 7712 , 771^) '[71,7171,7712 , 7713) 



d7l72 m2 

G3{n,7i7i, 7712,7713) = F3 (n, m^,m2,m3,p^ = -(mi + m2 - ma)^) 
— if mi, m2, ms) 

^ii^ ' '[71,7711,7712,7713) H ii ^ ' '[71,7711,7712,7713) 



87713 ^3 

(9) 

By inserting the above expansions in the four equations Eq.(^ and in the three equa- 
tions Eq.(0) for i = 1,2,3, we find only five algebraically independent equations for 
Ga{7i, 7711, 7712, 7713), a = 0,1,2,3 and their derivatives in respect to 7713, which can be 
written as 



(n-3)(3n-8) 

Go(n, mi, 7112, 7713) 



(71 - if 

+ (n — 3)mi(2mi — 7713 + m2)Gi(n, mi, m2, 7713) H — Tin, 77il)T(7i, mo) 

4m2m3 

(71 - 2)2 

+ (71- 3)m2(2m2 - 7713 + 7711)02(71, mi, m2, ms) + T(n, 77i\)T{n, ttiI) 

4mim3 

{71 - 2)2 

+ in - 3)m3(2m3 - mi - 7712)03(71, mi, m2, m3) Tin, 77i\)Ti7i, ml) , (10) 

4mim2 

7713 - 7711 - 7712 d 

7^ — Go(n, mi, ma, m3) = [7i - 3)Go(n, mi, m2, m3) 

2 C'm3 

+ 77ilOi{7i, mi, m2, m3) + 77il02{n, mi, m2, m3) + 7713(7711 + 7712)03(71, 7711,7712, 7713) (11) 
and, for i = 1, 2, 3, 

d 

(7713 - mi - 77i2)(77i3 - mi)(m3 - m2)T^ — Oi(7i,77ii, 7712,7713) = 

07713 

Pi,i(7i, mi, m2, m3)Gi(n, mi, m2, mg) + Pi^2(n, mi, m2, 7713)02(71, mi, m2, m3) 
+ Pi^3(n, mi, m2, m3)G3(n, mi, m2, m3) + (5i,i(n, mi, m2, m3)T(ra, m^)T(n, m^) 

+ Qj,2(^, mi,77l2, 77l3)T(7l, 7n\)T(7l, 77x1) + Qj,3(«, mi, m2, 77l3)T(7l, 77ll)T(7l, 77ll) , (12) 
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where T{n,m?) is defined in Eq.@, while the rational functions Pij{n,mi,m2,ms) and 
Qij{n, nil, ""^s) (which differ from the functions of the same name, but different argu- 
ments, appearing in 0]) are defined as 



and 



^1,1 

P2,2 
P2,3 
P3,l 
^3,2 



R 



3,3 



n,mi,m2,m3) 
n, mi,m2,m3) 



n — 3 



2 

n — 3 



■m2{m3 - m-i) - (ms - mi)(m3 - m-i) 



'ni2{m-i - nil] 



n, nil, m2, m^) = -{n - 3) (^^{m^ - nii - 7722) + ^m2{nii + 7712) j 



n, mi, 7712,7713) 
77,7711,7712,7713) 
n,nii,ni2,ni3) 



77 — 3 



2 

71 — 3 



mi (7113 — m2) 



mi(m3 — mi) — (m3 — mi)(m3 — m2) 



-(77 - 3) ^m3(m3 - mi - m2) + -mi (mi + m2)j 

71, mi, m2, m3) = -(71 - 3) — (m3 - m2) (7713 - -^mi - m2) 

m3 \ 2 J 

71, mi, m2, m3) = — (71 — 3)-^(m3 — mi) im-}, — mi m2^ 

m3 V 2 / 



7i,mi,m2,m3j 



m3 



-(m3 — mi)(m3 — m2){nii + m2 — 2m3 



+ 



71 — 3 



2m3 



(4(m3 — mi)(m3 — m2)(m3 — mi — m2) + mim2(mi + m2 — 2m3)) , (13) 



(5i,i(77,mi,m2,m3) 
Qi,2in,nii,ni2,nis) 
Qi,3{n,mi,m2,ni3) 
Q2,iin,mi,ni2,nis) 
Q2,2{n,nii,ni2,ni3) 
Q2,3i'f^,mi,ni2,nis) 



8m2m3 



(71 - 2f 



, 2 2"("^3 - "^l)("^3 - ^"^2) 



77 



],mini2 



■{ni2 - mi) 



(77 - 2f 



{n-2f ^ 

2"("^3 - mi) 



[n 



mimg 
2)2 



iniini2 



(mi - m2) 
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Q3An,mi,m2,m3,) = ^3 (ma - m2){m^ - ^mi) 



Q3,2(n, mi, m2, ms) = -^^^^^ — ^(ms - mi)(m3 - \m2) 

Arriiml 2 

(n - 2)^ 1 1 

Q3,3(n, mi, m„ m,) = - ^^^^^^^ (^3 - - 2"^2) • (14) 



The three Eq.([l^) form a system of equations involving only the master amplitudes 
Gi{n,mi,m2,m3) i = 1,2,3 and not Go{n, mi, 1712, 1713), while Eq.(|lOD can be regarded 
as an identity expressing Go{n, mi, m2, ma) in terms of the Gi{n, mi, m2, ms); that means 
that only three of the four master amplitudes are independent at the pseudothreshold and, 
as a consequence, Eq. (|ll|) is implied by the other four equations and is not an independent 
equation. 

As any system of k linear first order equations is equivalent to a single k-th order 
equation, we transform the system Eq.(^) into a single third order differential equation 
for a single unknown function, which we choose to be Go{n, mi, 7712,1713); further, we use 
Eq.([To|) and Eq. ([Tl|) for expressing the Gi{n,mi,m2,m-i) in terms of GQ{n,mi,m2,m-i) 
and its derivatives with respect to m^. Once the equation for GQ{n, mi,m2, m^), is solved, 
i.e. GQ{n,mi,m2,m3) is known explicitly, the Gi{n,mi,m2,m3) can be easily recovered, 
as we show later. 

Before continuing, we expand all the functions in n around n = 4; according to the 
notation of illl we write 



Ga{n,mi,m2,m3) = G (n) 

+G'^^\mi, m2, ma) + 0{n - 4) 



- — ^-T^G^^ ^)(mi,m2,m3) + ^)(mi,m2,m3) 
[n — 4)^ [n — 4) 



a = 0,1,2,3. (15) 



The function G{n) has been already introduced above; again from |1[], where the singular 
parts of Fa{n,m\,m\,m\,p'^) are given for arbitrary values of p^, we have 

G'[)"^^(mi,m2,m3) = -\{ml + ml + ml) , 

o 

G'[)"^^(mi,m2,m3) = - — (mi + ms - m-^f + — (m? + m2 + m3) 

32 lo 

- - [mi log(mi) + m2 log(m2) + ml log(m3) 

G'S"^\mi,m2,m3) = i , ^^"^^(mi, m2, ms) = --^ + ^ log(m- ) , z = l,2,3. (16) 

8 16 8 
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The third order differential equation for G'o'^^(mi, m2, ms) then reads 

+5(1 (mi, m2, m3)-^G'[)°^(mi, ms, mg) + go{mi,m2, m^)G^Q\mi,m2, m^) 
orris 

+g{mi,m2,m3) = , (17) 

were the Qairrii, m2, m^) a = 0, .., 3 are the following polynomials 

5(3(mi,m2,m3)=m3(m3 - mi){ms - ms)^^ 
g2{mi,m2,ms) = {m-i - mi){m-i - 7/12) 
gi{mi,m2, ^3)= — + 2(mi + m2)x'^ + (3m^ + 8mim2 + 3ml)x^ + 4mim2(mi + m2)x 



1/ 

5m3 + - (mi + m2j 



5(o(mi, m2, ^3)= - 6 ^m3 - ^("^i + "^2)^ (^rnl - ms{mi + m2) + ^mim2^ 



'1^ 



with X = m3— mi— m2, while the function g{mi, m2, tris), which is written in the Appendix, 
is a combination of similar polynomials and of first and second powers of the logarithms 
of the masses generated in the expansion in (n — 4) of the T{n,m?). In order to solve 
Eq.(|l^), we make the ansatz 

G'[,°^(mi,m2,m3) = a(mi, m2, m3)G[,°'^^(mi, m2, m3) (19) 

with the idea of fixing the otherwise undetermined function a(mi,m2,m3) by imposing 
that in the resulting differential equation for G^Q'^\mi, m2, m3) the coefficient of the term 
without derivatives vanishes, i.e. by imposing 

5f3(mi, m2, ms)^—^a{mi, m2, m3) + 5-2(^1, m2, m3)^^a(mi, m2, m3) 
d 

+ gi{mi,m2, ms)^ — a(mi, m2, m3) + go{mi, m2, m3)a(mi, m2, m3) = , (20) 
Cm3 

which of course is nothing but the homogeneous version of Eg . ([TTD . We further look for 
a{mi,m2, m^) in the form of a rational function in m3 and mi, m2; by trial and error we 
find that Eq.([T9|) can be given the explicit form 

r-(o)/ N 2m3-mi-m2 (0,1) 

Go{mi,m2,m3) = -— -^Gq ' {mi,m2,m3). (21) 

2(m3 — mi — m2)"' 



With the above substitution Eq. ([T7|) has taken the form of a second order differential 
equation in the first m3-derivative (the essential property of the transformation Eq.(pTD 
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is that Go°'^^(mi, m2, ms) does not appear anymore in the equation). By iteration of the 
method, the following chain of substitutions is obtained 



d 



dm- 



d 



dm-: 



G'[,°'^^(mi,m2,m3) 



-{21713 -mi- + ao + 



(2/77.3 — mi — 



GrK,m„m3) = ]-a2Gt'\mi^m2^m,) 



(22) 



where 



ao 
ai 

a2 



-{?>m\ + 2mim2 + ?>m\ 



6 



12 



{mi + 3m2)(3mi + 7712) (mi — 7/12) 



(2m3 - mi - m2)\Jm2,{m3 - mi - m2) 
(m3 — mi)2(m3 — m2)^[m3 — (mi + m2)m3 — |mim2]^ 



(23) 



G'o°'^''(mi, m2, m3), finally, satisfies the following first order differential equation 



d 



dm- 



'[mi,m2,m3) 



-- (m3(m3 - mi - m2)) 2 go + gi log(m3) + g2 log^("33 
+ (m3(m3 - mi - m2))"^ q_i om^^ + q^i^m^^ log(m3) 

+g-l,2"3j^ log^("33) + ^0,0 + 90,1 log(m3) + go,2 log^("33) 

+qi,om3 + qi,ims log(m3) + gi,2m3 log^(m3) 

+(l2,oml + q2,iml log{ms) + q2,2ml\og^ {m^) 

+<l3,oml + q3,iml log(m3) + ^3,2^3 log^(m3) (24) 

+q4,oml + q^^iml log(m3) + gs^om^ + gs.im^ log(m3) 

where the coefficients g^, qa^jj are functions of mi and m2 only and are reported in the 
Appendix. 

Eq. (|2^) is lengthy, but in fact it is a very simple differential equation, which can 
immediately be solved by quadrature. From Eq.(p2[) and Eq.(p^ is obvious that we need 
to integrate three times in m3 to get G^Q^ (mi,m2,m3). The integrations are elementary 
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and the integration constants can be fixed at = mi + 1712; indeed, we are considering 
the pseudothreshold value = — (mi + m2 — ms)^, which for that value of corresponds 
to = 0, where the master integrals of the sunrise graph are known in closed analytic 
form Once the integration constants are fixed, the explicit form of Go°^(mi, m2, ms) 
reads 



Gr(mi,m2,m3) = 
1 

8(m3 — mi — m2) 
1 

4(m3 — mi — m2) 



m'^(mi + 2m2)£(mi, m2, m^) + m\{2mi + m2)C{m2, mi, m^] 



(mi + m2) [■m\C{mi, m2, ms) + m\C{m2, mi, ms^ 



+ -{mi + 1712 + fnim2) I mi log I — I + m2 log I — 



ms 



mi 



ms 



m2 



+ -{ml - ml- ml) {C{mi,m2,m^) + £(m2,mi,ms)) 
o 

1 

~32 



m\ log^(m^) + m\ log^(m2) + m\ log^(m3) + {m\ + m^ — m\) log(m^) log(m2) 



+ {m\ — m\ + m\) log(m^) log(m3) + {—m\ + m\ + m\) log(m2) log(m3 
— mi(7mi + 2ms) log('^i) ~ m2(7m2 + 2ms) log(^2) 
+ (2m^ + 2mim2 + 2m2 — 5m3) log(mg 

--^(m^ + m2 + ml) - ^(-mim2 + mims + m2ms) 
12o u4 



(25) 



where 



£(mi,m2,ms) 



ms 



m2 



Lis 1 - — - Li2 ^ + lo^ 



mi 



m2 



ms 



mi + m2 



log ( — 

m2 



(26) 



The above result is in agreement with To check it we have simply to expand in n 
around n = 4 our overall coefficient C{n) ( Eq.(9) of [jl|) and the n-dependent overall 
coefficients of [^]). 

As already mentioned, once GQ{n,mi,m2,m^) is known in closed analytic form, its 
derivatives are also known, and the other amplitudes Gi{n, mi, m2, ms), G2{n, mi, m2, ms) 
and Gs{n,mi,m2,m3) can be easily obtained through simple relations with them. 

To obtain such relations, we differentiate GQ{n,mi,m2,m3) in its integral representa- 
tion with respect to mi and m2, as done already in Eq.([TT|) for ms, and we get after some 
algebra the following two additional relations 
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d 



■Go{n,mi,m2,m3] 



-2 



[n - 3)Go{n, mi, 1712, 7713) 



dmi ~ ' ms — mi — m2 

+ miim-i - m2)Gi{n, mi,m2, m^) + mlG2{n, mi,m2, m^) + mjGsin, mi, m2, ms) 



d 



dm' 



■Go{n,mi,m2,m3) 



m^ — mi — m2 



[n - 3)Go{n,mi,m2,m3) 



+ m\Gi{n, mi,m2, m^) + m2{m3 - mi)G2{n, mi,m2, m^) + mlG^in, mi, m2, m^) 



(27) 



The first of tlie Eq.(|I^), Eq.(|TUp and Eq.(P7|) form a system of equations, which al- 
lows to express master amplitudes Ga{n, mi, m2, m^) , a = 0, 1,2,3, as a function of first 
derivatives of GQ{n,mi,m2,m3) only. As GQ{n,mi,m2,m3) is already known, we can 
express Gi(n,mi,m2,ms),i = 1,2,3, as a function of Go{n,mi,m2,m3) and two of its 
first derivatives (we choose derivatives in respect to mi and m2 ). The solution reads 
(the expression for G2{n,mi,m2,m3) can be easily found from Gi{n,mi,m2,mz) by the 
substitution mi — > m2 , m2 — > mi ) 



Gi{n,mi,m2,m3) 
1 



Smim^ 
+ 



d d 
(mi - 3m3) - — Go{n, mi, ms, m^) + (ms + mg) - — Go{n, mi, ms, mg) 
ami om2 

2{n - 3)(mi + m2) - {n - A)m3^Go{n, mi, m2, m^) 



+ 



m^ — mi — m2 
n-2f ( 1 



2{n — 3) \m2 
G3{n,mi,m2,m3) 
1 



T{n,m2)T[n,m^) H T{n,mi)T[n,m^) T [n, mi)T [n, m2j 



mi 



mim2 



^ml 



d d 
(3mi - mg) - — G'o(n, mi, ma, mg) + (3m2 - mg) - — G'o(ri, mi, ma, mg) 
ami am2 



1 



mg — mi — m2 
in - 2f / 1 



(n — 3)(6mi + 6m2 — 7m3) — mg )G'o(?t-, "^i, m2, mg 



f28) 



2(?2 — 3) \m2 



T{n, ml)T{n, mg) H 7'(n, m^)T(n, mg) —T{n, ml)T{n, ml) 



mi 



mim2 



Expanding around = 4 we find that ^-'(mi, m2, mg) and ^-"(mi, ma, mg), (a = 
0, 1, 2, 3), already explicitly given in Eq. (pIB]) , automatically satisfy the Eq. (pSj) , while for 
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Gf^\'mi, 7712, ^3), {i = 1, 2, 3), we find the explicit expressions 
Gf\mi,m2,m3) = 

mi (mi + 2m2)£(mi, m2, m-^) — m\C{m2, mi, m-i) 



8(m3 — mi — m2)^ 
1 

4(m3 — mi — m2) 



(mi + m2)C{mi, m2, ma) + ^ ^mi log (^^^ + "^2 log 



/mg 
Vm2 



+ ^(/:(mi, 



m2,m3) + /:(m2,mi,m3)j + — [log^(m^) + log(m^) log(m^) 



log(mi) log(m3) — log(m2) log(m3) — 41og(mi) + 21og(m3 



1 



Gf(mi,m2,m3) 
1 



^(m3 — mi — m2)^ 
1 



m1C{mi, 1712, TTT's) + m2i2(m2, mi, m^) 



nil 



mi log — + m2 log — 



m3 



m2 



8(m3 — mi — m2) 

~(^C{mi, 1712, ms) + C{m2, mi, ms)) + ^ log^(m3) - log(mi) log(m2) 



+ log(mi) log(m3) + log(m2) log(m3) — 2 log(m 



(29) 



where C{mi, 1122, m^) is defined by Eg. (p6D . The function G2'\'rTii,m2, m^) is obtained 
by substituting mi m2 , m2 mi in the expression of ) (mi,m2,m3). 



4 The expansion at the pseudothreshold. 



As recalled in [Q, it is easy to find the coefficients of the expansion of a function in any 
variable once a system of differential equations in that variable is given. In this section 
we obtain the expansions in at the pseudothreshold of all the master integrals Eq.(0), 
Eq.(||). Putting expansions Eq.(|) into the system of equations Eq.(|) we get (comparing 
the lowest terms in expansion) from the first equation 

1 



n, mi, 1712,1713) 



[n - 3) G^in, mi,m2, m^) 



(30) 



(m3 — mi — m2)^ 

+ mlGi{n, mi, m2, m^) + mlG2{n, mi, m2, m^) + m\G'i{n,mi,m2,m'i) 
Substituting for the Gain, mi, m2, m^) their expressions expanded around n = 4, we have 

1 



H^^'^\n,mi,m2,m3) = G'^in) 



32 [n - 4) 



+ if(°'^)(n = 4,mi,m2,m3) + C(n-4) 



, (31) 
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with 



H^^'^\n = 4, nil, ""^2; ^3) 
1 



+ 



^(ma — nil — 1^2)^ 
1 

]{ni3 — mi — ni2Y 



8(777.3 — TTii — 7772) 



mi + r?72 

16(r?73 — 7771 — "^2) 



777^(7771 + 2r772)£(mi, 7772, "^3) + m\{2nii + m2)C{ni2, nii, nis) 
-2(7771 + "^2) {ni\C{nii, 7772, "^3) + mlC{ni2, mi, 7773)) 

, /777i7772\ 3, /J^T-A 3, /"^2 \ 

+ 777i 7772 ( 777i + 7772) iog ^ + 777^^ log + 7772 

\ 7773 y V?773/ \7773/ 

m\C{mi, 7772, ^^3) + 7772>C(7772 , 777i, 7773) — - [m\ + 777 i7?72 + 7772^ 



777i 



7773 



777i -777i + 7772 log - 777 2 777i + -"^2 Uog 



7772 



7773 



log(7773) ^ . 

16 ^ ^ 128 



(32) 



The other three equations form a system of hnear equations, which can be easily solved 
algebraically. The solution for H^^'^\n, mi, m2, m^) reads 



^(1,1), 



77, 7771, 7772, 7773) 



1 



mi{m3 — mi — 7772)" 



77 - 4 V 16 



— (7773 - 27771 - ^2) Gi{n, 7771, 7772, "^3) 



+ 777 (^3 - 7771 - 27772) ^2(77, 7771, "^2, 7773) 

16 

+ (27773 - nil- ^12) G^{n, mi, 7772, 7773) 
16 



(33) 



We do not put here the rest of the expression as it is not relevant for the subsequent 
discussion. Similar formulae are found for H^'^'^\n, 7771, 7772, 7773) and H^^'^\n, mi, 7772, 7773). 

As it is evident from the Eq.(^), to get terms ~ (77 — 4)" in the (77 — 4) expansion 
of H^^'^\n, mi, 7772, ""^3) we need terms ~ (77 — 4) of the Gi(77, 7771, 7772, 7773) , 7 = 1, 2, 3, (a 
similar problem was found in |^ for the expansion at the threshold). The terms ~ (77 — 4) 
are not yet known and could in principle be investigated with the method used in the 
previous section. The integrals obtained that way contain up to trilogarithmic functions 
(which are however expected to cancel out in the combination appearing in the r.h.s. of 
Eq.(^) - see below), so we decided to get H'^'^'^\n,mi,m2,m^) ,7 = 1,2,3, in a simpler 
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way. By using Eq.@ we get the following relations 

H^^'^\n, 7711,1712,1713) -- 

[71,7711,7712,7713) = 



dH^^'^\7l, 7711, 7712, 7713) ^ 2 m3r^i-^2 ^(0,2)^ 

+ 2- 



+ 2— -H'"'^' {71,7711,7712,7713) 

OTTLi mi 

dH^^^^\n, 7711, 7712, 7713) , o"^3-"^l-"^2^(o,2)/ 



^(3,1), 



71, TTli, 7712,7713) 



87712 7712 
[7l,77li,77l2,77l3 



71, 77li, 7712,7713) 



dni\ 



) ^mg-mi-ms (0,2)/ 



7713 



7712,7713) 



(34) 



As the pole terms in [n — 4) expansion of H^^'^\7i, mi, m2, ma) ,i = 1, 2, 3, are known to 
vanish from Eg . (13^) and with the result in Eq.(pl|) we deduce that H^°''^\n, 7711,7712, 7713) 
has no pole at n = 4; the Eq. (|3^ takes the form for n = 4 

dH^'^'^\7i = A, 7711,7712,7713) 



H^^^^\n = A,77ii, 7712,7713) = C^(ri) 



+2 



7713 — f^l — 



mi 



d7n\ 

^(0,2)^^ _ 4^ 



H^'^''^\7i = 4,7711,7712,7713) = C^in) 



dH^^'^\7l = 4, 77li,77l2, 77I3) 



+2 



ms — mi — m2 



7712 



^(0,2)^^ _ 



H^^^^){n = 4^^ mi, 7712,7713) = C^i^n) 



(0.1) = 4, mi , m2, ms) 
(9mq 



ma - mi - m2 ^(0,2) 



mg 



' (n = 4, mi, m2, ms) 



(35) 



Therefore to find H^'^'^\7i = 4, mi, m2, ^7^3), i = 1, 2, 3, it is enough to insert the derivatives 
in respect to the masses, easily obtained from Eq. (p2D , and to have the explicit form of 
(n = 4,7711,7712,7713). We evaluate it in the next section using dispersion relation 
method. 

There is no further obstacle in obtaining higher order terms in the expansion at 
the pseudothreshold, as can be verified with an explicit algebraic calculation, by solving 
iteratively the system of linear equations, obtaining the coefficients of any order in terms of 
the coefficients of lower orders. The reason is that the troublesome denominator 1/ (n — 4) 
appears only when expressing the first order coefficients in terms of the zeroth order 
coefficients; the denominators appearing at higher orders are l/{n — 5), l/(n — 6), ■ ■ ■, 
regular at n = 4. 
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The explicit expressions of the solutions are easily obtained with an algebraic program, 
but they are somewhat lengthy and we do not present them here for the sake of brevity. 
Let us also observe that for the practical purposes of solving the system of differential 
equations numerically only the values at pseudothresholds and first derivatives are needed. 

It is worth to note that the functions (^^(mi, m2, m^) ,i = 1, 2, 3 are not independent 
(the pole terms in Eq.(|5B[) have to cancel) and the relation between them, using Eq.(|TBp 
and Eq. (|29|) , reads 



+ 



+ 



+ 



+ 



Gf'\mi, 772,2, 272,3) (—8mi ~ 4r?7,2 + 42/23) 
G'^\mi, 1712, "2723) (—4mi — 8m2 + 4m3) 
G^\'mi, 2722, 2773) (—4mi — 4m2 + 82773) 

— 2?7i — 1712 + ^^3) 



1 r 
2 



1 r 



1 r 



+ log(?72^)(-m2 + 2773) + log(m2)(-mi + nis) - log(m3)(mi + 2272) 

log^(277^)(2mi + 7772 — 7773) 

+ log^(2772)(7r2i + 27772 — ""^s) + log^(r?73) (tT/i + 7^22 — 2m3) 
log(rr2^) log(2772)(277i + 1712) 

+ log(777^) log(r773)(r77i — 7^23) + log(r772) log(?723) (7^22 — 7^23) 



(36) 



This relation is of course fulfilled by the already known solutions Eq.(|29|). 



5 The value of H^^''^\n = 4, mi, 1712, ms) • 

Let us define the function F{p'^) through the equation 



Fo{n = 4, 777^, 7772, ^4, V^) = Fo{n = 4, 777^, 7773, 7773, 0) 

+ ^0(^^ = 4,7771,7772,7773,0) 



+ 



F{P' 



(37) 



The subtracted dispersion relation for F{p'^) reads 
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~. 2^ . 2^2 7 du 1 , 2 2 2 ^ 

Fip') = {P^f / ^, ^ 2^ - I"^^o n = 4, m , m , m^, u) , 3J 



(mi+m2+»n3)2 

where 



(^il-nia) 



1^^. .0 0-. 1 /■ „ Riu,h,m\) R(b,m'7,ml) , , 

- ImFo(n = 4, m?,m2,m2,n) = — j dh ^ ' ^ ^ ^' , (39) 



(mi+m2)^ 

with the usual definition of the function R 



y ^ z) = \J + y"^ + — 2xy — 2xz — 2yz . (40) 

The function H^^''^\n = A,mi,m2,m^) which we want to calculate in this section in 
the notation of Eq.(P7D,Eq. (P5|) reads 



H^^''^\n = 4,1711,1712,1713) = -F" {p^ = —{m^ — nil — 1712)'^) = 

I — — - \m.FQ{n = A,ml,ml,ml,u) , (41) 

J (m — [m^ — nil — m2) ) it 

(mi+m2+ma)2 

By using Eq.(|39|) and exchanging the order of integration we get 



i7(°'2)fn = 4 



,mi,ni2,nis) 



1 7 R{b,nilnil) J RMmj) ^^^^ 



16 J b J u {u — {ni^ — nil — "^12)^)^ 

(mi+m2)^ (Vb+m^y 

The integration in u is elementary and gives 

H('''\n = 4,mi,m2,m,) = ^ J dh R^^^^l^D H{h) , (43) 



(mi+m2)^ 

where 
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lZ{h)(m^ — mi — 1112] 



2 



2{ml - hf 2(6 + ml) Abrnj 



+ 



(m3 — mi — m2)^ (ms — mi — m2)^ lZ^{h) 



+ log 



h \ —h + m 



2 



ml y (m3 — mi — m2)^ 



with 



and 



2 1 6 + m? 1 

H 1 , (44) 

(m3 — mi — m2)'' (m^ — mi — 1712^ TZ'^{b) TZ'^{h) ' 



7^(6) = i?(6, (mi +m2 -m3)^m^) (45) 

^ _b + ml- (mi + m2 - m3)2 - TZ{b) ^^^^ 
2 m3 



To perform the last integration, it is convenient to differentiate with respect to one of 
the masses Eg. (^31) in its integral form, to perform the integration in h of the derivative 
so obtained (the integration in h being now elementary) and then to integrate again in 
the mass. As a final result we find 



rT(o,2)( . \ ^ \ R'^imi,m2,-m3 

' '[n = A,mi,m2,m3) 



32 [ 2 (m3 — mi — m2)^ 

2 mi log 



mi 
ma 



(m3 — mi — m2)^ 



2 m2 log 



(ml + mg — mim2 + mim3^ 
(ml + mg — mim2 + m2m,^ 



(m3 — mi — m2)^ 

+7 ^ 2mi mn — m.mn — m^m^ — — m, — m^ 

(m3 - mi - m2)'^ V ^ ^ ^ ^ ^7 (7723 - mi - m2)'5 ^ ^ 

4 X2 mim2m3 



(m3 - mi - m2)^ R{{mi + m2)2, (mi - m3)2, (m2 - m3)2) j 
where the integrals Xj i = 0, 1, 2 are defined in the following way 



(47) 



db 



(mi+m2)^ 



log (1/5) 



+ 



log 



(Vb 

Vm3 



i? (6, (mi — m3)2, (m2 — fn^y) R{b, mf, m^) 
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= (,^ks('^)los{^)-U,(l-^)-U,{l-^) , (48) 



d6 



(»ni+m2)2 



(mi — m2){mi + m2 — 2m3) 
R (b, (mi - m3)2, (ma - mg)^; 



log(|/5) 



(m^ — 7712 



2 log 



-C2 + log ( — ) loj 

\m2. 



17121713 



{nil + "^2)' 



2Lb I 

m2 



L12 ( 1 

ms 



i?(6, mf, 1172 
Li2 



ma/ 



(49) 



and 



I9. 



db 



R ((mi + ma)^, (mi 



"^3)^ im2 - m^y) 



{mi+m.2)'^ 



5C2 + log 



- (mi + ma) 
'tl (l + t?)\ 



i? (6, (mi - mzY, (m2 - m^Y) 



'^2 



log 



.2 .2^ 

'-1 '-2 



/2 

^2 



- 2U2 
+ U2 

where 



t2(l-tl 
. t2-l 



2 Lio 



+ Lis 



l-t2 



- log 

- 2 Li2 

Li2f 



l-t2 
1 + tl t2 
1 +t2 



log 



- 2 Li 



(1-^2)^ 
tl + t2 



l+t2 



^2 ^ 

1 + 



1 + 



(50) 



y/mi + 7172-1173 - Jml 



i/mi + m2 - ms + i/mi 



and to 



/mi 



mi + ^m2 



(51) 



The last expression was found assuming mi > m2 > ma, but the analytic continuation 
to other regions is straightforward. 



6 Summary 

In this paper we have presented the expansion of the 2-loop sunrise selfmass master am- 
plitudes at the pseudothreshold = —(mi + m2 — ms)^. The other pseudothresholds can 
be easily found by the permutation of the masses. We define the expansion in Eq.(^; the 
values of the amplitudes at the pseudothreshold are given in Eq. (|15]) , Eq. ([T6| ) , Eq. (pSj) and 
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Eq.(^). The first order terms in the pseudothreshold expansion of the master amplitudes 
at {n = 4) are presented in Eg. (pT]) , Eq.(|32D and Eq.(^), while the second order term 
of FqIti = 4, m^, m^, mg,^^ = —(mi + m2 — ^3)^) is presented in Eg. (^Tj) . The higher 
order terms, which are not given exphcitly here, can be easily found by solving recur- 
sively, at each order, a system of four algebraic linear equations. The expansion at the 
pseudothreshold cannot be simply deduced from the known expansion at the threshold 
and vice versa, even if at first sight they seem to be connected by the change of sign 
ms — >■ —ma. In fact the analytic properties of the amplitudes are different at the two 
points: at the pseudothreshold the sunrise amplitudes are regular, so that the solution of 
the Eq.(^) can be expanded as a single power series, while at the threshold the sunrise 
amplitudes develop a branch point and its expansion is indeed the sum of two series, . 
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7 Appendix. 

In this appendix are given the definitions of the functions used in the text. 
Function g{mi,m2,rn^) used in Eq.(|T7p 

5((mi,m2,m3) = Tq + n log(m^) + log(m^) + rg log(m^) + ri^i log^(m?) 

+■^1,2 log(mi) log(m2) + ri,3 \og{m\) log(m3) + r2,2 log^(m2) 
+r2,3 log(m^) log(m^) + r3,3 log^(m^) , (52) 

where 

11 . 147, , 4 ,485 2 211 485 3 

ro = -—X - ^("^1 + ^2)x - + -^^1^11^2 + —m^jx 

,81 o 839 2 839 , 81 2 
lb o4 o4 lo 

— ( — mi + 5mim2 H m^m2 + bminin H rriojx 

16 2 16 

7 

— — mim2(m'^ + 2m\m2 + 2mim2 + Tn\) 
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ri 

^2 
^^3 



mi 



ml 



1.9 9 5 5 3 

o^^^ + o ("^1 + "^2)3;^ + (TT^rnj + -mim2 + 777"^2)3; + T7:"^i"^2("^i + "^'2) 
2 o io 4 ib ib 

1 9 5 5 9 3 

15 5 71, ,4 ,121 2 67 121 2x 3 

— X H mi + m2 X + m, H mim2 H m, x 

4 8 ^ ^ ^16 ^ 4 16 

+ (3m? + -—rrri\m2 + — — mim, + 3m2)x^ 
lb lb 

+ (— m^ + —m\m2 + —m\m\ + —mim\ + — mg)^ + — mim2(mi + m2)^ 



n,2 
^1,3 

^^2,2 
^^2,3 
^3,3 



3 

16 



m\ ^ms - -^("^1 + "^2)^ [rn\ - m^{mi + m2) + ^mim2^ 
^mlml (^ms - i(mi + m2)^ 

— ml (m^ (mi + m2)'j (Sx"^ + 3(mi + m2)x + m2(mi + m2) 
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— m2 ma - -(mi + m2) ) [m^- m^{mi + m2) + -mim2 



— -m2 ^ms — 2 ("^1 "I" "^2)^ (3x^ + 3(mi + m2)x + mi (mi + m2)^ 
(^{ml + mg)^^ + 3(m'^ + m\m2 + mim\ + m2)x + mim2(mi + m2)^) 



and X = ms — mi — m2. 

The definition of tlie coefficients used in Eg . (|2^) : 



[7 1 1 

go = mim2(mi + m2) — ("^-i + mim2 + m\) - —m\ log(mi) - —m\ log(m2) 

'-48 lb lb 

+ log(mi) - m2 log(m2))^ 

9b 

qi = m\m\{mi + m2)^ " « ("^1 + "^2) + 777"^! log(mi) + —m2 log(m2) 

L 8 24 24 



?2 
9-1,0 



■^mlml{mi + m2)^ 



mim2 
1 



192 



(5mi + 29mlm2 + 104m^m2 + 29miml + 5m2) 



+ — m^(m^ + 5m2) log(m^) + —nT'l^bml + log(m2) 

+ — mim2(mi log(m^) - m2 log(m2))^ 
48 



(53) 
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7 1 
m\ml{mi + 777,2) "^("^i + "^2) + (7771 log (7771) + 7772 lo, 

5 5 21 4 67 3 2 67 2 3 21 

+777^(^777? + ^777^7772 + ^mi777^ + ^777^) log(777?) 
+777^(^777^ + ^"^?^2 + ^^l"^2 + ^^2) log(777^) 

rn\rn\{rrii + 7^72) log(m?) log(m2) 

16 



777 i 7772 
1 



+ ^777 i 7772 ( 777i + 7772) (log(777i) + log(7772)) 

o 

-777 ^7772 (7771 + 7772) 

65 4 1031 o 857 2 2 1031 , 65 . 

777i H 777?7772 H 777i777, H 777i7779 H 777, 

64 ^ 192 ^ 96 ^ ^ 192 ^ 64 ^ 

+ ^"^1^2 + Y6"^2) l0g(777i) - ml(—'^l + ^"^l"i2 + Yq'^I) l0g(^2) 

— — (7771+7772)^ 777^ log^ (777^) + 7r72 log^ (7773) H 777^7772 log(777^) log(7772) 

32 L J 8 

A 137 Q 25 o <-) 137 o 4 

777^ + ^^777^7772 + —Tn-^m2 + -^mim^ + 7772 

-777i(^777i + ^77717772 + ^7773) log(777i) - 7772(^777i + ^777 i7772 + ^7773) log 

— - (^3777^ + 67/7^7772 + 10777i7r72 + 617111712 + 37r72j 

,rl93 , 1775 

-(7771 + "72) 
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777^ 



"61^^ + 192 



193 2 

-77717772 + -gj"^2 



^777^ log(777^) + ^777^ log(777^) - ^777^ log^(777i) - ^7772 log 

63 



777, 



(7771 + 7772) - ^"717772 - y7772 + ^777i log(777i) + ^7772 log(7772) 

9 , X / O 0\ 



-{mi + 7772) (777i + 7772 
o 

251 2 241 251 2 

-^rn, + —77717772 + —7772 



b4 'Z4: b4 

+ -777ilog(777i) + -7772 log(7772) " ^^"^1 Wl^^l) " ^^"^2 log^("^2) 
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^3,1 = + 39mim2 + -p"^2 - ^"^i log("^i) - ^"^2 log("^2) 

93,2 = -^{ml + m^) 

37, , 79, , 11 15 

94,0 = -^("^1 +^2) g4,i = — -(mi + m2) ^5,0 = -7-7 95,1 = ^ (54) 
32 4 lb 2 
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